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Motivated by the recent experimental observations on clustering of motor proteins on microtubule 
filament, we study an open system of two parallel totally asymmetric simple exclusion processes un¬ 
der asymmetric coupling conditions, which incorporates the mutual interaction with the surrounding 
environment through Langmuir Kinetics in both the lanes. In the modified Langmuir Kinetics, the 
attachment and detachment rates depends on the configuration of nearest neighboring sites. We 
analyse the model within the framework of continuum mean-field theory and the phase diagrams 
along with density proHles are obtained using boundary layer analysis. The effect of mutual inter¬ 
actions on the phase diagram for two different situations of attachment and detachment (LK) rates 
is discussed. Under the symmetric LK dynamics, the topological structure of the phase diagram 
remains similar to the one in without mutual interaction; while for the antisymmetric case, after 
a certain critical value of attractive/repulsive mutual attraction, signihcant changes are found in 
the qualitative nature of phase diagram. Moreover, it is shown that the type of mutual interaction 
affects the dynamic properties of motor proteins. The theoretical hndings are examined by extensive 
Monte-Carlo simulations. 

PACS numbers: 05.60.-k, 87.16.Nn, 64.60.-i, 05.70.Ln 

I. INTRODUCTION 

Many natural systems exhibit complex behavior un¬ 
der stationary state when either driven by some external 
field or self driven. Such driven diffusive systems reveal 
very rich nonequilibrium phenomena in physics, chem¬ 
istry and biology such as kinetics of bio-polymerization 
[l| , dynamics of motor proteins in biol^ical cells [3, 3 ; 
gel electrophoresis vehicular traffic [sj, and modeling 
of ant-trails Q. In biological systems, molecular mo¬ 
tors are motor proteins that consume chemical energy 
and move along polymer filaments of the cytoskeleton, 
which act as macromolecular tracks. Their collective dy¬ 
namics plays a major role in various intracellular pro¬ 
cesses/functions such as cellular trafficking, protein syn¬ 
thesis, cell division etc iai- 

In order to analyze the collective properties of interact¬ 
ing molecular motors, totally asymmetrically simple ex¬ 
clusion process (TASEP) model is found to be a paradig¬ 
matic model to study motion of motor proteins in the last 
decade [§-[ll|. TASEP model comprises of single species 
of particles performing biased hopping with uniform rate 
in a preferred direction along a ID lattice. The parti¬ 
cles obey certain preassigned rules under hard-core ex¬ 
clusion principle, due to which a lattice site cannot have 
more than one particle. In an open TASEP connected 
to boundary reservoirs, the particles remain conserved in 
the bulk. The conservation of particle number is violated 
by allowing particle adsorption (desorption) to (from) the 
bulk with fixed rates [Langmuir kinetics (LK)] similar to 
the attachment-detachment of motors occur form the fil¬ 
ament [m . Single-channel TASEP coupled with LK has 
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been well studied [isl - fTsIl and the phase behavior is con¬ 
siderably different from those known in reference models 
of TASEP without LK Q. The competing dynamics of 
particle conservation (TASEP) and particle nonconserva¬ 
tion (LK) results into distinguishing characteristics such 
as localization of shocks and phase coexistence ©[a. 
The bulk and surface transitions in a single TASEP with 
LK are examined via boundary layer analysis [l^ . 

Experimental studies on collective transport of kinesin- 
1 motor protein along microtubule suggests that the 
individual motors interact with each other via short 
range interaction and this attraction is weakly attrac¬ 
tive [IMl- The simulation results of the motor attach¬ 
ment/detachment dynamics in the presence of mutual in¬ 
teractions are found compatible with in-vitro experimen¬ 
tal observations [II- Moreover, this mutual interaction 
manifest as a propensity of the motors to cluster during 
their transport on microtubule. As a result, intermolecu- 
lar interactions affect various chemical transitions inside 
the cell and play an important role in the collective dy- 
namics of molecular motors. The results in refs, [l^, [2j 
further strengthen the fact that mutual interactions be¬ 
tween motors plays a significant role in determining the 
spatial organization of the motors on filaments. To quali¬ 
tatively understand the collective transport of weakly in¬ 
teractive molecular motors through short-range interac¬ 
tions on single filament, two different approaches namely 
mean field and modihed cluster mean held approach is 
emp loyed on TASEP model with [2^ and without LK 
1^ . respectively. Recently, the effect of mutual interac¬ 
tions is examined in a one dimensional TASEP model by 
modifying the attcahment/detachment rates while keep¬ 
ing the hopping rate unaffected [1^. This is in contrast 
to the Katz-Lebowits-Spohn (KLS) model [2^ which in¬ 
clude both the correlations due to mutual interactions 
and hard core exclusion principle. 
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It is evident that more realistic description of many 
transport processes such as motor protein, vehicular traf¬ 
fic etc. can be mimicked via multi-channel TASEPs. 
Moreover, kinesin (a motor-protein) can move along par¬ 
allel protofilaments [2^ forming a multi-lane system. 
Due to its importance, the two-channel system is studied 
extensively without LK under different coupling environ¬ 
ments [26l - l28l| . Recently, the competition between bulk 
and boundary dynamics in steady state is examined by 
Gupta and Dhiman [2^ in a two-channel TASEP model 
with Langmuir kinetics in both the lanes under biased 
lane-changing rule. But the important aspect of mutual 
interaction is ignored in this study. 

In this paper, we explore the consequences of mutu¬ 
ally interacted LK on the stationary density profiles as 
well as on phase diagrams of a two-channel open TASEP 
under the fully asymmetric coupling environment. The 
proposed framework not only provide natural means to 
consider the interaction between the motors with the sur¬ 
rounding environment at micro level but also be useful 
to understand some aspects of the collective behavior of 
motor proteins that emerges at macro level in a coupled 
system. We also investigate the role of symmetry of inter¬ 
action via LK rates and show that it significantly affect 
the system dynamic. The paper is organized as follows: 
In the next section, we firstly present the theoretical de¬ 
scription of the two-lane TASEP model under the biased 
coupling environment. To incorporate the mutual inter¬ 
actions with the environment, the modified version of LK 
is proposed followed by the mean field approximation. In 
sec III, the method to obtain stationary density profiles is 
discussed. The phase diagram and density profiles under 
two different cases of modified LK rates are investigated 
in sec. IV. Finally, in sec. V, we provide our conclusions. 


II. TWO-LANE MODEL AND 
HYDRODYNAMIC MEAN-FIELD 
APPROXIMATION 

We define a two-lane open system of two parallel one¬ 
dimensional lattice channels, each with N sites, denoted 
by A and B. The system consists of indistinguishable 
particles distributed under hard-core exclusion principle 
and hop unidirectionally to the right (See Fig. [T]). The 
system configurations are characterized in terms of occu¬ 
pation numbers (i = 1, 2,3,. N] j = A, B), each of 

which is either zero (vacant site) or one (occupied site). 
In an infinitesimal time interval dt, a lattice site (i,j); 

i = 1,2,3,.V; j = A, B is randomly chosen and the 

paticle hoping takes place according to the following dy¬ 
namical rules. 

1. If z = 1 (entrance of lane j), particles are injected 
into the lattice with a rate a provided tij = 0. 
If Ti j = 1 and r 2 j- = 0, then the particle moves 
from site (l,j) to site (2,j) with a unit rate. If 
Ti j = 1 and r 2 j- = 1, no hopping occurs from 


site (1, j). Neither lane changing nor attachment- 
detachment takes place at z = 1. 

2. If z = V (exit of lane j), a particle is removed 
from the lattice with a rate /3 when r/vj = 1- 
Neither lane changing nor attachment-detachment 
takes place at z = V. 

3. For 1 < z < V (bulk of lane j), if = 0, a 

particle can attach to the site with a rate (at¬ 
tachment rate). When ti^a = 1, then the particle at 
the site (z, A) firstly tries to detach itself from the 
system with a rate (detachment rate) and if it 
fails, then it moves forward to site (z -|- 1, A) with 
unit rate provided = 0. If a particle at site 

(z, A) fails to detach from the system followed by 
an unsuccessful attempt of moving forward, then 
it will try to shift to lane B with a rate w; pro¬ 
vided Ti^B = 0. Here the parameter i depends on 
the configuration of the neighboring sites and is a 
non-negative integer which represents the strength 
of the modifying factors of attachment/detachment 
rates. The dynamics in the bulk of lane B are simi¬ 
lar, with the only exception being that particles are 
forbidden to shift from lane B to lane A. Here, no 
lane-changing is allowed from lane B to lane A. 

Here the mutual interactions are incorporated by modi¬ 
fying the Langmuir Kinetics in terms of detachment (at¬ 
tachment) rate from (to) the background. The modified 
LK dynamics at site (z, j) depend on the configuration of 
two nearest neighboring sites according to the following 
rules: 

• If both the neighboring sites are vacant i.e. j = 
0 and Ti+ij = 0, then the parameter £ = 0 resulting 
the attachment (detachment) rate as Wa (wd)- 

• If either of the neighboring sites are vacant i.e. 
Ti-ij = I and Ti+ij = 0; or Tt-ij = 0 and 
Ti+ij = 1, £ takes the value 1 and the attachment 
(detachment) rate becomes jWa (Swd)- 

• If none of the neighboring sites are vacant i.e. 
Ti_i J = I and Ti_|_i J = 1, attachment (detachment) 
rate adapted as 'j^Wa {S'^Wd) due to £ = 2. 

For 7 = (5 = I, the proposed model reduces to a two- 
lane asymmetrically coupled TASEP with LK in both the 
lanes [^. The rules for modified LK rates in both the 
lanes imparts the model its generality over the existing 
two-lane coupled TASEP models with LK [2^, [s^. Note 
that our model is suitable to study a number of two-lane 
transport processes such as vehicular traffic and motor 
proteins because normally a vehicle or a molecular motor 
does not change its lane unless hindered by another one 
preceding it. This is completely in accordance with the 
lane-changing rules defined by us. 

The master equations for the temporal evolution of 
occupancy of particles at each bulk site (1 < z < Y) in 
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both the lanes (j = A, B) is given by 

d’^Tij'^ /i \ /i \ 

B,j)^ ^2+1, 

“t“Wa<C(l Ti_i^j)(l T2 j')(1 

+7u;a<(l - Ti7)[(l - ri^ij)ri-ij + (1 - ri-ij)Ti^ij]> 

+7^Wa<ri_lj(l - Tij)Ti+ij> 
-UJd<Ti,j{l - Ti_ij)(l - Ti+ij)> 
j'((1 “t“ (1 

( 1 ) 


where < • • • > denotes the statistical average and last 
term on right-hand side takes a negative (positive) sign 
for lane-A (B). At boundaries {i = l,iV), the occupancy 
of particles evolve according to 


d<Ti,j> 

dt 

d<TNj> 

dt 


= a<(l-n.i)>-<rij(l 
= <TAr-lj(l — TN,j)> — 


-T2,j)>, 

( 3 <TN,j>. 


( 2 ) 

(3) 


In the above system of equations correlations are neglect¬ 
ing by using mean-field approximation consist of the fol¬ 
lowing approximation. 


<Tjjri+ij> Ri <Tij><Ti+i7>. (4) 

For large system size N ^ oo and small lattice spacing 
e = L/N -A 0 with finite iVe, one can derive the con¬ 
tinuum limit by coarse-graining a discrete lattice with 
quasi continuous variable x = iL/N and rescaling the 
time as t' = tL/N. For simplicity, the length of the lat¬ 
tices is fixed as L = 1 which restricts the variable x in 
the range [0,1]. The time rescaling is useful to under¬ 
stand the engagement between particle conserving and 
non-conserving dynamics as the system attains station¬ 
ary state locally due to conservative dynamics only and 
the non-conserving processes occur at a comparatively 
lower rate than particle conserving processes. To study 
the competition between boundary and bulk dynamics. 


we define the reduced the attachment, detachment and 
lane-changing rates inversely proportional to the system 
size (T^. [i^. 


Xla = ujaN, Xld = uJdN, XI = ujN. (5) 


Note that the parameters Xla-, Xld and XI remain fixed as 
N -A oo. 

In the continuum limit average density is replaced by a 
continuous variable by substituting <Tij> = pij G [0,1] 
and retain the terms up to 0{N~‘^) in Taylor’s series 
expansion (for large system i.e. N » 1) to obtain 


_ I 1 ^Pi,j 
P^.j±i - P^,o ± ^ 


1 d'^pij 

2Af2 dx'^ 



( 6 ) 


Without loss of generality, we drop the subscript i as 
the system is free form any kind of spatial inhomogene¬ 
ity. The continuum mean field equations describing the 
steady state average densities (pA and ps) in the bulk is 
given as 


A 


PA 

d 

“h —— 

Pb_ 

dx 


+ PAi^ ~ pa) 
m+PBa-PB)\ 


= S, 


(7) 


where 

■na(l - Pa)(1 + Pa{S - 1))^ - ^dPA{^ + Pa{i - 1))^' 
-^Pa^^-Pb) 

f^a(l - Pb)(1 + Pb{5 - 1))^ - XtdPBi^ + Pb{i - 1))^ 

L +xip\{i-pB) 

Here, the term S originates from the combina¬ 
tion of particle lane changing transitions and modified 
adsorption-desorption kinetics and is responsible for the 
loss of particle conservation in the bulk of each lane. The 
components of J represent the currents in the particle 
conservation situation in lane-A and B, respectively. The 
important noteworthy aspect arising due to the biased 
lane-changing rule is that the coupling term acts as a 
sink for lane-A and a source for lane-B [2^. Addition¬ 
ally, the particle-hole symmetry does not hold for the 
model with mutual interaction. 
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III. STEADY-STATE SOLUTION 


In this section, we analyze the coupled system obtained 
form mean-field description under steady-state and dis¬ 
cuss the effects of important parameters on density pro¬ 
files and phase diagrams. The system dZD in the steady- 
state reduces to 

2 ~ ~ Pa){^ + Pa{5 — 1 ))^ 

+ PA (7 - 1))^ - ^^Pa(1 - Pb) =0^^ 

2 ^^2 ■*" ^‘^P^ ~ ~ Pb)(1 + Pb{S — 1))^ 

—^dPsi^ + Pb{i — 1))^ + ^^Pa(1 ~ Pb) = 0- 


The boundary conditions for coupled nonlinear system 
(0) are p^(0) = ps(0) = a, pa(1) = Pb( 1) = 1-/3. In 
the continuum limit e 0 , the coupled system becomes 
overdetermined as the second order nonlinear differen¬ 
tial equations reduce to first order differential equations 
keeping all the four boundary conditions intact. Here, 
the leading order terms in the above system is analogous 
to the dissipative term of viscous Burgers’ equation in 
stationary state and play a similar role as performed by 
the vanishing viscosity term in the solution of Burgers’ 
equation. 

To understand the steady-state behavior of the system 
o, we employ boundary layer analysis on the contin¬ 
uum mean-field equations. Recently, this approach has 
been quite successful in solving the hydrodynamic equa¬ 
tion in the thermodynamic limit and explaining the com¬ 
plete rich phase diagrams of single-channel TASEP with 
LK [I^ as well as two-channel TASEP with LK in fully 
and partially asymmetric coupling conditions [ 2 ^[3l| , re¬ 
spectively. The density prohles are constructed by com¬ 
puting the bulk part of the solution (outer solution) and 
boundary layer solution (inner solution) separately and 
then to match the both solutions suitably to get the 
global solution. The outer solution is the major part of 
the density profile and can be obtained in the thermody¬ 
namic limit [L » 1) due to the negligible contribution 
of the regularizing terms. We use a suitable numerical 
scheme to get approximate outer solution of the con¬ 
tinuum mean-field equations. The following numerical 
scheme for lane is used to find the outer solution of 
continuum mean-field equations in both the lanes. 


n-l-l _ 

Pi.j 


' 2 Aa ;2 

At' 


2Ax 


{Pi+l,j ‘^Pi,3 3 " Pi-1,j) 

[( 2 p((^--i)(pr+i.,-pr-i,,)] 


+ At'[B„(l-p((^.)(l + a,(<5-l))2 
- + p,.,(7 - 1))^ T f^(p- a)^(1 - pIb)] 

+ 0{At',Ax^). 


Here, the last term takes negative (positive) sign for 
lane-A (B), respectively. The stationary state density 
profiles have been obtained by capturing the solution of 


above system in the limit n >> 1 , which ensures the oc¬ 
currence of steady state. This solution is also refereed 
as bulk solution. The solution describing the bound¬ 
ary layers or shocks can be obtained by ignoring the 
non-conservative terms in the steady-state system. This 
can be achieved by using a transformation x = , in 

the e-neighborhood of another position variable Xd repre¬ 
senting the position of boundary layer in the continuum 
mean field equations. This leads to the elimination of the 
source and sink terms in the hydrodynamic equations in 
regions of width of 0 (e) and is well justified because par¬ 
ticle non-conserving dynamics are irrelevant in regions of 
width of 0(e). For more details of this approach we refer 
to [H,[11,133. 


IV. PHASE DIAGRAMS AND DENSITY 
PROFILES 

We restrict our investigation for a special choice Ha = 
VLd to simplify the analysis and construct the phase dia¬ 
grams. Importantly, the biased coupling of lanes makes 
the steady-state dynamics nontrivial in contrast to a sym¬ 
metrically coupled two-lane TASEP with LK which is 
similar to two independent single-lane TASEP with LK 
dll. In order to investigate the effect of mutually inter¬ 
active LK dynamics on the steady-state properties, the 
whole discussion is divided into two parts based on the 
symmetric or antisymmetric LK rates. In the symmet¬ 
ric case, attachment and detachment rates are modified 
in a similar fashion while in the antisymmetric case, if 
attachment rate is enhanced than the detachment rate 
is reduced by the same amount and vice-versa. The re¬ 
sults obtained from continuum mean-field equations are 
validated with Monte-Carlo simulations utilizing random 
sequential update. In the simulations, a lattice of size 
L = 1000 is used to minimize boundary effects and the 
simulations are carried out for 10^*^ time steps. To ensure 
the occurrence of steady-state, the first 5% steps are dis¬ 
carded and the densities in both the lanes are computed 
by taking time averages over an interval of IQL. The 
phase boundaries are calculated within an estimated er¬ 
ror of less than 1 %. 


V. CASE 1: MUTUAL INTERACTION WITH 
SYMMETRIC LK RATES 

In this section, we construct the density profiles and 
derive phase diagram for B = 1 and = 0.2 under the 
mutually interactive symmetric LK dynamics. To inves¬ 
tigate the effect of mutual interactions with symmetric 
modified LK, we choose 7 = ^ = 1 -|- 0 i.e. both the at¬ 
tachment and detachment rates are symmetrically modi¬ 
fied with a multiple of 1 - 1-0 for each occupied neighboring 
sites with 0 > — 1. Here, 0 represents the strength of the 
mutual interaction under symmetric LK rates. It is to be 
noted that the positive (negative) values of 0 increases 
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FIG. 2. (Color online) Phase diagrams for = 1. (a) 0 = 0 
i.e. no mutual interaction, (b) 9 = 0.5, (c) 9 = 2.0, (d) 
9 = —0.5, (e) 9 = 0 . In (a)-(d): fid = Ha = 0.2 and in (e): 
Hd = Ha = 0.5. Here, LD, HD and S denotes Low density. 
High density and Shock phase, respectively. 


(decreases) the LK dynamics provided the neighboring 
sites are occupied. 

For 0 = 0, the present case will convert in to a simple 
asymmetrically coupled two-lane TASEP with the LK 
system without mutual interaction whose steady-state 
behavior has been recently investigated [ 2 ^ and found 
to be quite complex as compared to the symmetric cou¬ 
pling. It has been reported that a small asymmetry in 
lane changing rates can produce significant changes in 
the phase diagram without mutually interactive LK . 
Figure 2(a) shows the phase diagram in terms of bulk 
transitions with 0 = 0 , which has been adapted from 
ref. [ 2 ^ and is provided here for the sake of compari¬ 
son. Since the surface transitions are responsible only 
for the change in the slope of the boundary layer, we 
for simplicity, focus on the simpler representation of the 
phase diagram showing only the bulk transitions in the 
phase plane. Moreover, this simpler structure depicting 
the qualitative nature of phases is helpful in analyzing 
the effect of mutual interactive LK on the phase dia¬ 
gram. Broadly classifying, there exists six steady-state 


distinct phases viz. (LD,LD), (LD,S), (LD,HD), (S,S), 
(S,HD), and (HD,HD) as shown in fig. [^a). Due to the 
asymmetric coupling between the lanes, lane A acts as 
a homogeneous bulk reservoir for lane B. This creates 
an imbalance between adsorption and desorption rates 
in both the lane leading to a higher effective detachment 
(attachment) rate for lane A (B). As a result, major re¬ 
gion of the phase diagram is filled with LD (HD) phase 
for lane A(B). This combined effect can be seen in the 
phase diagram of individual lanes A and B which are sim¬ 
ilar to the one channel TASEP with LK for and 

Hq > Od, respectively. 

To investigate the effect of mutual interaction with the 
background, the phase diagrams for nonzero 0 are shown 
in figs. [2Kb)- (d). It is clear from the figures that the (HD, 
HD) phase disappears quickly from the phase diagram 
while the (LD, LD) phase shrinks slowly on increasing 
0. This is due to the fact that for higher 0 there will be 
more detachments (attachments) in HD (LD) phase due 
to the higher probability of occupied neighboring sites. 
For 0 = 0.5, the number of steady-state phases reduces 
to five, consisting of all the phases except (HD, HD) as 
compared to the one with 0 = 0. For 0 = 2, phase 
(LD, LD) becomes confined to a very small region near 
to the boundary a = 0 (fig. HKc))- The reverse phe¬ 
nomenon happens when we start decreasing 0 from zero. 
The phases (LD, LD) and (HD, HD) expand while the 
other phases shrinks. Note that there does not emerge 
any new phase for negative values of 0 . 

Under the symmetric LK rates, the effect of 0 on the 
phase diagram is found analogous to the one with Hq. To 
compare the effect of 0 with Ha, the phase diagram for 
higher value of Ha = 0.5 for 0 = 0 is drawn in fig. He). 
The topology of the phase diagram for increasing 0 is 
similar to the one corresponding to increasing Ha with 
the only exception that (LD, LD) phase still present in a 
narrow region for higher 0. The main difference between 
varying 0 and Ha is that the impact of 0 is much stronger 
in the (HD, HD) phase as compared to Ha. 

Fig. |3| shows the density profiles of all five different 
phases for 0 = 2.0 from continuum mean-field equations 
and their validation using Monte Carlo simulations. In 
general, there is a good agreement between the simu¬ 
lation and theoretical results except near the boundary 
layer due to the finite size effect. It is to be noted that 
PA < Pb because of the fully asymmetric coupling as 
particles of lane-B cannot move to lane-A. The effect of 
0 on (LD, LD) phase is shown in fig.IlK^-)- From the vari¬ 
ation of density profile with respect to 0 , it is clear that 
the phase (LD, LD) disappear with an increase in 0 and 
converts in to mixed type of phase in lane B for larger 
values of 0. The mixed type of phase corresponds to the 
two-phase coexistence of LD with HD. It is to be noted 
that the LD-HD phase is not characterized by a discon¬ 
tinuity or a shock localized in the bulk of lane B. This 
is due to the fact that the LD and HD profiles are not 
straight lines for 0 > 0 which is in contract to the case 
of 0 = 0. Moreover, the density profile in both the lanes, 
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(a) (a) 




FIG. 3. (Color online) Density profiles in symmetric case with 

6 = 2.0 and Q,d — fla = 0.2 (a) (LD, LD) for ot = 0.03 and 

7 = 0.3; (LD, S) for a = 0.2 and 7 = 0.3; and (LD, HD) for 
a = 0.8 and 7 = 0.3, (b) (S, S) for a = 0.2 and 7 = 0.65; and 
(S, HD) for a = 0.8 and 7 = 0.65. The solid (dashed) lines in 
red (blue) color are the continuum mean-field density profiles 
of lane A (B). The curves marked with squares (circles) are 
the result of Monte Carlo simulations for lane A (B). 


regardless of the phase, approach to Langmuir isotherm 
as 6 increases. In the limit of 0 —>■ 00 , the Langmuir 
Kinetics dominates the system dynamics completely and 
the density in both the lanes becomes Langmuir density 
i.e. PA = Pb = 1/2. Further, the effect of 9 on the (LD, 
HD) phase is examined in fig.SKb). The existence of (LD, 
HD) phase for larger 9 is shown. In contrast to the case 
of 0 = 0, density becomes uniform throughout the lattice 
in both the lane for higher values of 9 and approaches to 
Langmuir density. 

A shock is a discontinuity in bulk connected by a low 
and high density part. Across the boundary layer in the 
bulk, i.e. a shock, the current in the system remains 
constant [l^. The effect of 9 on the shock dynamics 
in the (S, S) phase is explored in Fig. [SKa). The first 
thing one can notice is that 9 has a negative effect on the 
shock hight. The indistinguishable feature of the shocks 
in our system is that shock in lane A continuously moves 
towards right boundary on increasing 9 while the shock 
in lane B shows a different trend known as jumping effect 


FIG. 4. (Color online) Effect of 9 on density profiles in both 
the lanes for 0 = 0, 2, 5,10 (a) a = 0.03 and 7 = 0.3, (b) 
a = 0.8 and 7 = 0.3. The continuum mean-field and Monte 
Carlo simulation results are represented, respectively, by solid 
(dashed) lines in red (blue) color and squares (circles) for lane 
A(B). 

. As shown in fig. [5Kb) shock in lane-B moves to the 
right first till 0 < 1.5 and changes its direction for further 
increase in 9 and starts moving towards left boundary for 
9 > 1.5. Moreover, for comparatively larger value of 9 
shock in both the lanes disappears and bulk density of 
the whole lattice approaches to Langmuir density. 


VI. CASE 2: MUTUAL INTERACTION WITH 
ANTISYMMETRIC LK RATES 

So far, we have investigated the effect of modified LK 
rates when both attachment and detachment rates are 
enhanced or reduced simultaneously by equal amount 
while the hopping rate remains unchanged. Recently, in 
an in-vitro experiment, the presence of attractive inter¬ 
actions has been seen in kinesins-I, which remain longer 
attached to the microtubule in the presence of neighbor¬ 
ing motors and leads to formation of motor clusters la¬ 
in this section, the proposed model is used to understand 
the coordination mechanism under the existence of at¬ 
tractive interactions. The effect of mutual interaction is 
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(a) 




FIG. 5. (Color online) Effect of 6 on shocks in (S, S) phase 
for O' = 0.2 and 7 = 0.8 (a) 61 = 0, 0.5, 1.0, 1.5, 2.0, 5.0. Inset 
shows verification via Monte Carlo simulation results for N = 
1000, (b) Location of shock (xa) vs. 9 in both the lanes. 




FIG. 6 . (Color online) Phase diagrams in antisymmetric case 
with lid = fla = 0 . 2 , and = 1 for different values of cj>- (a) 
(f) = 0.5, and (b) 4> = —0.5. 


analyzed in antisymmetric manner where the attachment 
(detachment) rate is enhanced (reduced) and vice versa. 

To begin with we set j = 1 + (j) and S = 1 — (j). Here, (j) 
is a constant having the range [-1,1] as any value outside 
this interval leads to negative rate which is unrealistic. 
Positive values of (j) represents the attractive interaction 
while negative values signify the repulsive interaction of 
the LK dynamics. The effect of (f) on the density pro¬ 
files are analyzed and phase diagrams are constructed 


using the methods discussed in previous section. Firstly, 
we investigate the role of attractive interactions on the 
phase diagram by increasing the value of (j). For (j) = 0, 
the system has been well studied and have six distinct 
phases (see fig. [2][a))[2^. The increase in the value of cj) 
firstly retains the topological structure of the phase dia¬ 
gram with six distinct phases while the phase boundaries 
are slightly shifted. Due to enhanced attachment rates 
via mutual interactions, the region containing (LD,LD) 
phase shrinks, while the region of (HD, HD) phase ex¬ 
pands with increasing (j). The significant changes in the 
topological structure of the phase diagram are observed 
for (j) = 0.5 as shown in Fig. UKa). Note that there does 
not emerge any new phase and the phase diagram consists 
of five phase having the same characteristics as for (/> = 0. 
The (HD, HD) phase captures majority of the region in 
the phase diagram while the (LD, HD) phase completely 
disappears. In fig. [7l[a), the density profiles for three dif¬ 
ferent phases from continuum mean held along with the 
Monte Carlo simulation results are shown corresponding 
to the phase diagram in Fig. [6l[a) with D = 0.2. 

The transition from (LD, HD) to (HD, HD) with 
increasing (j) can be understood as follows (see hg. [8|). 
Due to the increased number of particles attaching from 
the background in both the lanes, the density in both 
the lanes increase with (/) which leads to abundance of 
particles in the system. As a result, very less number of 
particles shift from lane A to lane B for higher (j). An 
noteworthy aspect of the density prohle is the origination 
of a mixed type of profile for intermediate range of (j>, 
where some part of the bulk is in low density while the 
remaining part have high density. So, it is reasonable to 
conclude that LK dominates the lane changing dynamics 
in the presence of higher attractive interactions with the 
background. 


Now, we discuss the important and distinguishing fea¬ 
tures of the phase diagram for negative 4>. Though the 
repulsive interactions has not been studied in literature 
for motion of molecular motor, yet the effect of nega¬ 
tive values of (j) on the phase diagram has been explored 
for the sake of completeness as well as the results may 
be useful for explaining many other nonequlibrium phe¬ 
nomena. As contrast to the case of attractive interac¬ 
tions, (LD, LD) phase expands and covers major portion 
of the phase diagram for smaller values of (j). This can be 
understood on the similar lines as discussed for the case 
of (HD, HD) phase for positive (p. In fig. EKb), the phase 
diagram for (p = —0.5 are provided. There are only four 
different phases in the phase diagram which were also 
seen for the case of no mutual interactions i.e. (p = 0. It 
is to be noted that the topological structure of the phase 
diagram is quite different as compared to any other phase 
diagram discussed in this paper. Fig. [7] shows the den¬ 
sity profile of different phases computed form continuum 
mean field theory and compared with Monte Carlo simu¬ 
lation results. Although some discrepancies in the some 
part of bulk of Lane B is reported in Fig. [71 yet, there is 


























(a) 



FIG. 7. (Color online) Density profiles in antisymmetric case 
(a) (LD, LD) for a = 0.05, /3 = 0.4; (S, HD) for a = 0.3, fi = 
0.4 and (S, S) for a = 0.05,/I = 0.9 with (j) = 0.5, (b) (S, S) 
for a = 0.05, /3 = 0.9 and (S, HD) for a = 0.4, /3 = 0.99 with 
0= -0.5. 


a good agreement in general between the analytical and 
simulation results for antisymmetric LK dynamics. The 
discrepancies in Lane A occurs only near the boundary 
layer in the bulk i.e. shock and is due to finite size effect. 
On the other hand, whenever there is a shock in lane A, 
the discrepancies occur in the left part of Lane B form the 
position of the shock in Lane A. Moreover, the simula¬ 
tion results overshoot (undershoot) the analytical results 
for positive (negative) (p in the portion of discrepancies 
in Lane B. There are two reasons for such discrepancies, 
the finite lattice size and the mean field approximation. 
Moreover, the phase transition from (LD, HD) to (LD, 
LD) is shown in Fig. [5] Further, the effect of (p on the 
domain wall in (S, S) phase is analyzed in Fig. [HI In 
contrast to the symmetric case, shocks in both the lanes 
moves towards left boundary with increasing (p. 

The above analysis reveals that the mutual interacting 
LK dynamics has an important effect on the steady-state 
behavior of the two-channel system. Due to the recent 
findings on clustering of motor proteins on microtubules, 
this effect has been modeled by modifying the LK dy¬ 
namics based on the nearest neighbor interactions. 



FIG. 8. (Color online) Phase transition from (LD, LD) to 
(HD, HD) for a = 0.5, /3 = 0.5 with respect to an increase in 
(p. 



FIG. 9. (Color online) Effect of p on shock position in both 
the lanes for a = 0.05, (3 = 0.9. 


VII. CONCLUSION 

In this work, we study a simple two-lane TASEP model 
for the transport of molecular motors, which not only in¬ 
teract with each other but also with cellular environment 
through adsorption/desorption dynamics under a biased 
lane-changing rule. Based on the recent results on mo¬ 
tor protein on microtubules, the mutual interactions are 
incorporated by the means of modified LK dynamics de¬ 
pending on the configuration of nearest two neighboring 
sites. To explore the consequences of mutual interactions 
on the steady-state properties of the system, mean-field 
equations in the continuum limit are computed and the 
coupled system is analyzed using a singular perturbation 
technique. In general, the results of continuum mean- 
field equations agree reasonably well with Monte-Carlo 
simulations. 

The phase diagrams are obtained for two different cases 
of mutual interactions. Under the symmetric LK dy¬ 
namics in which both attachment and detachment rates 
increase or decrease simultaneously, the topology of the 
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phase diagram remains qualitatively similar to the one 
obtained in the case of without mutual interaction. The 
only changes in the structure of the phase diagram are 
the gradual shifting of the phase boundaries and the 
shrinkage/expansion of various phases. We have also 
analyzed the effect of mutual interactions on motion of 
shocks in the bulk of both lanes, their positions, and 
heights. For antisymmetric LK dynamics, it is observed 
that the topology of the phase diagram changes signif¬ 
icantly with an increase in attractive/repulsive mutual 
interaction. The effect of varying mutual interactions on 
the shocks in both the lanes are different as compared to 
the one with symmetric LK rates. 

The proposed work is an attempt to provide a nat¬ 
ural means to qualitatively understand the steady-state 


properties of a unidirectional transport on a two-channel 
lattice with mutually interactive LK under fully asym¬ 
metric coupling environments. The present study gives 
some insight not only in understanding complex dynam¬ 
ics of clustering of motor proteins but also towards en¬ 
hancement of one’s insight about various non-equilibrium 
systems present in nature. 
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